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La pr6sente invention conceme un procede et un recepteur 
pemiettant I'identif ication autodidacte des paramfetres caract§ristiques 
d'un canal de propagation k partir notamment des statistlques d'ordre 
2 du signal re9U par un capteur. 

Les paramfetres identifies sont par exemple le retard et 
I'attenuation des trajets du canal. 

L'Invention trouve son application notamment dans le domaine 
des radiocommunications mobiles ou des liaisons ionospheriques en 
HF. 

Canal speculaire 

• Les radiocommunications mobiles ou Ionospheriques HF sont 
sujettes au ph6nom6ne de trajets multiples. En zone urbaine, le signal 
6m\s est reflechi et diffract§ sur des obstacles fixes ou mobiles 
presents dans I'environnement Dans le cas des transmissions HF, les 
reflexions se font sur les differentes couches de la ionosphere. Dans 
le cas de canaux de propagation, pouvant §tre qualifies de speculaire; 
i.e. la transmission s'effectue selon un nombre limits de trajets discrets 
ou ponctuels caracterises par un retard et une attenuation complexe. 
En considerant un temps d'observation compatible avec la duree de 
stationnarite du canal, une reception multicapteurs et une propagation 
speculaire, le canal multi-capteurs a pour expression : 

c(t) = X)a4.J(t-'n:) W 

k=l 

oO k est I'indice d'un trajet, ak le vecteur dont les composantes sont les 
attenuations complexes du trajet k pour les differehts canaux, xk. le 
retard associ6 k ce kieme trajet et d le nombre de trajets dans un 
canal. 
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En outre, si chaque trajel k est Incident sur ie r#seau suivant-un cone 
de diffusion spatiale r^duH, {'expression (1) prend la forme sulvante 

d 

c(*) = Ii;/3*a(Wt-^) (2) 
ik=l 

ou a(ek) est Ie vecteur d'arrivee pour un angle associe au ki^me trajet 
et pk I'attenuatlon complexe du trajet. 

Dans Ie premier cas (Eq. 1), Ie canal n'est pas parametr^ par la 
direction d'amvee des trajets, chaque trajet de propagation est 
parametre par un temps d'aniv6e x et un « vecteur d'antennes » a. 
Ainsi, la calibration n'est plus n^cessaire et les algorithmes ne sont 
plus llmites par la dispersion spatiale ou les trajets coherents. 

Dans Ie deuxieme cas (Eq. 2), les trajets sont parametres par 
leurs directions d'arriv6e, ce qui suppose que la vari6t6 de Tantenne 
est connue et entratne done en gendral la mise en oeuvre d'une 
calibration pour estimer les valours de Sk- 
Identification autodidacte 

Dans les systemes actifs ou entrames, les parametres du canal 
sont calcules au cours d'une phase d'apprentissage durant laquelle 
r^metteurtransmet une sequence connue du recepteur. 

Si Ie canal de propagation fluctue au cours du temps, 
notamment a cause des mouvements des stations mobiles ou des 
couches ionospheriques, la sequence doit etre emise periodiquement 
af in de mettre ^ jour la yaleur des parametres. 

Si de tels systemes sont efficaces, remission r6guliere d'une 
sequence d'apprentissage peut entratner une reduction importante du 
debit efficace. 

Par exemple, dans la norme STANAG 4285 pour les 
transmissions HF ionospheriques cooperatives la moiti6 des syniboles 
transmis sont des symboles d'apprentissage. 
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L'art anterieur decrit aussi differentes methodes et systemes 
autodidactes, ou systemes aveugles, dans lesquels, restimation des 
parametres est effectuee h partir des statistiques du signal regu, sans 
connaissance a priori de la sequence d*apprentissage. 
Techniques d'ordre deux 

Les techniques propos^es utilisent, par exemple, uniquement. 
les statistiques d'ordre deux (matrice de covariance spatio-temporelle) 
du signal regu. Par rapport a des algorithmes aux ordres superieurs, 
les algorithmes a I'ordre deux ont de meilleures proprietes de 
convergence, i.e. pour un nombre de symboles donne, la variance des 
estimateurs a Tordre deux est en general inferieure a celle des 
estimateurs aux ordres superieurs. En autre, elles presentent moins 
de problemes locaux des problemes d'optimisation que les techniques 
aux ordres superieurs. 

De nombreuses methodes et algorithmes ont ete developpes, 
tels que ceux d^crits dans le document ayant pour titre "Multichannel 
Blind Identification : from subspace to maximum likelihood methods." 
de LTong et SPerreau ; Proceedings of the IEEE, 86(10) : 1951- 
1967, October 1998. Uun des Inconvenients de ces algorithmes est 
qu'ils nScessitent la connaissance de la longueur du canal de 
transmission global et que la plupart ne sont pas du tout robustes a 
una erreur d'estimation de cette longueur. Or, lorsque la transmission 
est a bande limitee, la longueur du canal n'est definie 
qu'approximativement et ces algorithmes ne peuvent etre utilises. 
Etat de Tart en identification parametrique autodidacte 

Les principaux travaux d'identification autodidacte des 
parametres caracteristiques de la propagation sont des algorithmes 
decouples, i.e. estimant independamment les temps d'arriv^es et les 
vecteurs d'antennes ou les directions d*arrivee. On peut citer les 
algorithmes de goniometrie pour Testimation des directions d'arrivee, 
les algorithmes de separation de sources (coherentes) pour 
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I'estimation des vecteurs d'antennes. II existe aussi d© rrombreux 
travaux concemant restimatlon des temps d'arrivSe. 

Les methodes d'estlmatlon conjointes conduisent a amiliprer la 
precision et la resolution des estimateurs permettant par exempio, 
I'estimation des parametres m§me lorsque les retards (ou les angles) 

sont tres proches. 

Les principaux travaux d'estimation conjointe autodjdacte des 
parametres (e.x) sont presentes, par exemple, dans les references 
suivantes : 

• "Identification spatio-temporelle de canaux de propagation a trajets 
multiples", de J.Gouffraud, PhD thesis, Ecole Nomfiate Sup^rleure de 
Cachan, 1997, ou 

• "Improved blind channel identification using a parametric approach" 
de M.C.Vanden/een et A.Paulraj, IEEE Communications Letters, pages 
226-228, August 1998. 

L'idee est que les entires utilises pour I'estimaWon liutodidacte 
de la reponse impulsionnelle peuvent etre minimises directement en 
fonction des angles et des retards, par des criteres de type sous- 
espace, tels que decrits dans la reference "Subspace methods for the 
blind identication of multichannel FIR filters" de E. i\/loulines, P. 
Duhamel, J-F.Cardoso et S. Mayrargue, IEEE Trans.on signal 
Processing. 43(2) : 516-525, February 1995. Ces algorithmes 
necessitent la connaissance a priori du filtre d'emission/reception et 
une calibration de I'antenne. 

D'autre part, un exemple de methods d'estimation conjointe 
autodidacte des parametres (a,i;) est par exemple decrit dans la 
reference "Methods for blind equalization and resolution of overlapping 
echoes of unknown shape, de A.Swindlehurst et J.Gunther, IEEE 
Trans.on Signal Processing, 47(5) : 1245-1254, May 1999. Les 
auteurs travailient dans le domaine frequentiel et proposent un 
algorithme iteratif de type IQML pour approximer la solution du 
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maximum de vralsemblance ainsi qu'un algorithme d'initialisatlon 
expllcite base sur I'algorlthme ESPRIT. Ces algorithmes ne 
niqessitent pas la connaissance a priori du filtre d'emlssion, mais 
presentent I'lnconv^nlent de n6cesslter une transfomiation de Fourier 
avant traiterhent. 

L'invention concerne un proc4d6 autodldacle ou partlellement 
autodidacte permettant de determiner au moins les param§tres spatlo- 
temporels caracteristiques d'un canal de propagation dans un systeme 
comportant au moins un capteur de reception de signaux y(t) 
caracteris§ en ce qu'il comporte au moins une etape ou la stmcture de 
type spiculaire du canal est prise en compte et une etape de 
determination conjolnte des parametres spatio-temporel, tels que les 
vecteurs d'antenne (a) et/ou les vecteurs temporels (x) k partir des 
statistlques d'ordre 2 des signaux re9us sur un ou plusieurs capteurs. 

Selon un mode de mise en oeuvre, il comporte une etape ou le 
signal regu y(t) est sur6chantillonn§, par exemple avec une periode 
d'echantillonnage correspondant a T/p, ou T est la p§riode symbole. 

Selon un autre mode de realisation, le proc6d6 comporte une 
etape de reception sur au moins deux capteurs et une etape de 
surechantillonnage. 

Le proced§ selon invention trouve son application, par 
exemple, dans le cadre des liaisons normalis6es a des fins 
d'^galisation, de localisation et de filtrage spatial ou dans le cas de la 
sun/eillance du spectre d'un canal de propagation, a des fins de 
localisation k partir d'une ou de plusieurs stations HF. 



L'invention presente notamment les avantages sulvants : 
• il n'est pas n^cessaire de transmettre une sequence 
d'apprentissage pouvant §tre la cause d^une diminution du debit 
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efficace de transmission, ou de connattre la longueur du cJanal de 

transmission globale, ^ 

. elle s'avere robuste a la surestimation de I'ordre du canal de 

tra^nsmisslon. 

D'autres avantages et caracteristiques de I'lnvention 
apparaltront h la lecture de la description se referant a la figure unique 
donnee a tltre illustratif et nullement limitatif qui represente le schema 
d'un r6cepteur de radio-communication selon I'invention. 

Avant de d6crire les etapes du proc6d6 selon I'invention 
quelques rappels. hypotheses et definitions sont donnees. 

Hypotheses et module 

Le modele de signal mono-capteur 

Le signal emis est de la fomr»e 

a;(t) = - JT)s(i) (3) 

ou g(t) est le filtre d'emission/reception et {s(l)> sont les symboles 
d'information qui ont et§ emis a la periode-symbole T. 

Le systfeme suppose que la frequence porteuse et la periode 
symbole ont 6t§ estlmees de fa9on pr6alable. Le signal re9U en bande 
de base s'ecrit 



y(t) = ^ h(t - lT)s{l) + w(t) 



(4) 



ou {8(1)} sont les symboles d'infomiation qui ont 6te 6mis a la periode- 
symbole T, h(t) est la r^ponse impulsionnelle du canal de 
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transmission, et w(t) est un bmit additif. Le signal regu ecliantillonne a 
la p6riode-symbole suit le modele discret : 

y(fc) = y{kT) = X) - + «'(*) 

(ez 

ou h(l) = h(IT) (en supposant que la phase d'echantillonnage soit egale 
a 0). 

le systeme multi-voies 

La figure 1 represente un systeme de reception multi-voies 
comportant un nombre q de capteurs de reception et un moyen 
pemiettant de surechantillonner le signal s(k) conrespondant au trajet 
d'indicek. 

Sans sortir du cadre de I'invention, il est bien entendu que le 
proc6d6 conceme aussi des recepteurs comportant plusieurs capteurs 
sans moyen permettant de surechantillonner ou encore un seul 
capteur et un moyen de sur§chantlllonner le signal regu. 

Dans I'exemple de realisation donne ci-apr6s, le disposltif 
consld6r§ est un systdme multi-voies combinant une diversit6 spatiale 
et une diversit6 temporelle. I.e. le signal est regu sur un r§seau de q 
capteurs Ci ayant chacun une reponse h(i)(t) et il est ensuite 
surechantillonn4 par rapport k la periode symbole T, au rythme T/p 
comme represente sur la figure 1. Le sur^chantillonnage est realise a 
un moyen d'un disposltif adapte. La reponse impulsionnelle globale du 
canal de transmission et !e vecteur des observations multi-capteurs 
sont d^finis respectivement par les vecteurs colonnes qxl suivants : 
h(t) = [h(i)(t) h(q,(t)r (6) 

y(t) = [y (i)(t)...., y(q)(t)f (7) 

L'indice inf^rieur indique le numero du capteur et I'exposant T 
represente la transposee d'une matrice. Les caract6res gras indiquent 
une dimension spatiale. Les echantillons obtenus par echantillonnage 
au rythme T/p sont vectorises en rassemblant les signaux 
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echantillonn§s sur la m§me phase d'echantillonnage a rint^rieur d'une 
periode-symbole. Plus pr6cisement, a chacune des p differentes 
phases d'echantillonnage au rythme symboie est assocl6e une voie. 
creant ainsi sur chaque capteur un canal multi-voles a p sorties. On 
obtient ainsi un canal multi-voies a p.q sorties. La repons^ 
impulsionnelle globale du canal de transmission et le vecteur des 
observations multi-vpies sont d§finis par les vecteurs colonnes p.q x1 
suivants : 

h(k) = [h<^UkT)... h^k-nf (8) 

m = [y^'^ (kT)... (9) 
L'indice sup§rieur Indique la phase d'echantillonnage, 

h « (kT) = h(kT +(i - 1)T/p). y ^ (kT) = y(kT + (i - 1)T/p). 
Hypotheses 

Le precede selon I'invention utilise certaines hypotheses 
notamment les suivantes. 
Source de symboles 

H t - Le processus { s(k) } est par example un bruit blanc complexe 
qui peut aussi etre circulaire au second-ordre de variance unlt6, 

E[s(k)l = 0. E[s(k)sG)'^] = 8(k-j). Els(k)s(j)l = 0 (10) 
Bruit 

H 2 - Le bruit w(k) est un processus stationnaire blanc spatialement et 
temporellement, Gaussian, circulaire, de variance , 
E[w(k)w(k)"l = o2|pq C*"") 

avec Ipq la matrice identite 

Flltre d'emlssion/r6ception 

H 3 - Le filtre d'6mission/r6ception g(t) a un support f ini, g(t) = 0 

pour t [0, LgT[. oil LgT correspond a la longueur du filtre et T la 

periods d'echantillonnage. 

H 4 - Le filtre d'emission/reception g(t) <-> G(f) est approximativement 
a bande limit6e, de bande Bg 
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G(f)^0 V/f^[-B,.BjJ (12) . 

et G(f) ne s'annule pas dans une certaine bande 

fi^mf^O V/ . .. _ (13) 

Canal de propagation 

H 5 - Le signal est bande etroite pour le reseau de capteurs, I.e. on 
suppose que le temps de propagation du signal d'un capteur au 
suivant est tres Inferieur k I'inverse de la bande du signal. 
H 6-Le canal est constant sur I'lnten/alle d'observation. Le decalage 
Doppler et les r^sldus de porteuse sont supposes negligeables. 
H7' L'ecart tennporel entre le premier et le dernier trajet est fini, et on 
en connaTt une borne superieure : At^ . 

d etant le nombre de trajets et et ak le retard et le vecteur 
d'antennes de dimension q x 1 associes au kieme trajet la reponse du 
canal s'exprime de la fa9on suivante. 

d 

En §chantillonnant k T/p, 

oo i d 

hz{z)= J2 HkT/p)z-i' = J2^k9zirk;z) (15) 

°" j;x(r;^)^X)5(n--r)z-" (16) 

" n P 

est la transfomnee en z de filtre d'emission/reception decale de r 

et echantillonne a T/p. 
Par ailleurs, h(z) correspond k la reponse en z au rythme symbole du 
systeme combin6, 
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k=-oo 

Le signal re9U suit le module discret Equivalent mono-source multi- 
sorties : 

y(k) = [h(z)]s(l<) + w(k) (18) 
ii est interessant de noter que 

jhW(2) = Efcl-oo*»^*n*)'Z"* sont les composantes polyphases du 

canal surechantillonne. 

D'aprfes les hypotheses H3 et H7, la rEponse impulsionnelle du 
canal de transmission global muW-capteurs est a longueur finie. Done, 
les r6ponses discretes h(z) et hT/p(z) sont de degr6 fini. Nous notons L 
= deg(h) i.e. au moins une des voles a son degr6 egal a L, les autres 
pouvant etre de degre inferieur, 

h(z) = h(0) + h(1 )z ^ + ... + h(L)z-^ 
De meme, nous notons Lj/p = deg(hT/p). On a pL < Lr/p < p(L + 1)-1. 
Afin de simplifier la presentation, nous supposons que 

Lt/p + 1 =P(I-+1)- 

Par ailleurs, les algorithmes a I'ordre deux sont bases sur la 
condition de diversite du canal de transmission global suivante ; 
H 8 Le vecteur polynomial h(z) est irr^ductible, i.e. ces composantes 
n'ont pas de racines communes. 
Quelques notations et outils . 

Soient {b^^^ (z),..., b^Hz)) r polynomes scalaires, 

b(z) = [b (z) b^'^(z)f est un vecteur polynomial (un vecteur 

dont les coordonnees sont des polynomes) de dimension , r. On 
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appelle ordre ou degre du vecteur, deg(b(z)), le degre maximum des 
polynSmes 

[b^^\z), b<'\z)}. Si U est le degre de b{z). alors 



b(z) = X;b(fc)z 

fc=0 



-fc 



(20) 



Oil b(k) = [b<'^ {k) ..... b^'^ (k)f . Le vecteur associe a b(z). i.e. 
contenant I'ensemble de ses coefficients est designe a I'aide d'une 
fl^che, 



b = (b(0)^...b(L6y) 



(21) 



est de dimension r(Lb + 1) x 1. Done, en particulier h con^spond 
au vecteur de dimension 
q.p(L + 1) x1 associe a h(z) et hT/p{z) 
i.e. contenant I'ensemble de leurs coefficients, 

h = lh(0)^...h(L)T ^^^^ 

= (hx(o)^...hx(Lx)T : - '^^^ . 

De plus, soit K un Entler correspondant au nombre d'observations 
multidimensionnelles, K(b) est la matrice de Sylvester (I.e. Toeplitz 
par bloc) de dimension r(K + 1) x (K + U + 1) 



associ^e au polyndme b(z), et d^ftnie par 



rA'(b) = 



b(0) ... MLb) 0 
0 b(0) b(i^fc) 

0 ... 0 b(0) 



0 

b(L6) 
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Enfin, soit B(z) une matrice polynomiale de dimension r x D et de 

degree Lb , 

B(z) = [bi(z),..., bD(z)l 
ou bi(z) est un vecteur polynomial de dimension r. Nous notons 

B = (bi,...,b£,l 

la matrice associie de dimension r(U + 1) x D et Dk(B) la matrice de 
Sylvester par bloc de dimension r(K + 1) x (K + + 1)D definie par 

Vk(B) = [TKihi), .... Tk(hD)] (25) 

Xmin(A) designe la plus petite valeur propre de la matrice A et vpmin (A) 
le vecteur propre associ§ k Xfnin(A). 
Precede selon I'invention 

Le precede selon I'invention. repose notamment sur la prise en 
compte de la structure de canal de type sp6culaire pour detemnlner les 
param&tres specifiques d'un canal de propagation ou de transmission 
d'un signal notamment dans un systeme de communication tel que 
d6crit a la figure 1. Les paramfetres d6termin6s conjointement sent par 
exemple le vecteur temporel (t) des vecteurs retards (xk) et le vecteur 
d'antenne (a) des vecteurs (br ). 

Uhe premiere variante de mise en ceuvre du precede, ou 
"precede a priori", est particulierement bien adaptee lorsque le f litre 
d'emission/reception est connu. 

Une seconde variante de mise en ceuvre du precede ou 
"precede sans a priori" correspond mieux au cas ou les 
caracteristiques du filtre d'emission/reception ne sent pas connus. 

Precede avec a priori 

Le proced6 suppose que le filtre d'emission g(t) est connu. 
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Les differentes etapes du precede ou algoiithme avec a priori selon 
rinvention utilisent le fait que Tequation (15) s'ecrit sous forme 
"vectorielle" de la fa9on suivante : 

h = G.,.^(x)a (26) 

h est le vecteur de dimension q.p.(L + 1) contenant Tensemble des 
coefficients de la reponse hy/pCz), a est le vecteqr de dimension q,d 
contenant Tensemble des vecteurs d'antennes a = [af . . aj)^ 
et GL.d(r) est une matrice de dimension q.p(L + 1 ) x q.d contenant les 
versions retard^es echantillonnees k T/p du f iltre d'emlssion : 



( ^(O-n) ... ^(O-rj) \ 

^(f-n) ... i7(f-Trf) 



01, (27) 



p 

Iq est la matrice identity de dimension q. Nous notons x le vecteur de 
dimension d contenant I'ensemble des retards x = [ti,...,t^ ] 

Les differentes etapes du precede sont par exemple les suivantes : . 

1 . Choisir la longueur du canal global de propagation tenant compte 
de remission et de la reception, L > + ^^max ©t la valeur 
du nombre d'observatlons K>L avec At^ la plus grande valeur 
possible du retard relatif entre deux trajets. 

2. Calculer la matrice Qi(R) du critere sous-espace par exemple 
de la maniere suivante : 

- Calculer Testimee empirique de la matrice de covariance : 



fc=0 



(28) 



j?K(ifc) = WfyKvivf..-y{KT-TivYf (29) 

= (y(Of ...y(/0%(fc+i) 
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et N-K le nombre d'observations ind^pendantes. 

- Constmire la matrice 6^ en prenant par exemple les vecteurs 
propres assocl^s aux pq(K + 1) - (K + L + 1) plus petltes valeurs 
propres de R- peut aussi etre obtenue a partir de la 
decomposition en elements propres de la matrice des donnees ^ 

y ^ [f(^K) . . . y(iv)].) Puis former la matrice de projection = G^Gf . 

- Former la matrice du critere sous-espace Qf,(R) = Vi{Ui)Vt(J^t) ■ 

3. Estimer le nombre de trajets d par exemple en mettant en oeuvre la 
methode decrite. par exemple dans la reference "Detection of 
signals by infomiation theoric criteria" de M. Wax et T. Kailath, 
IEEE Trans.on AcoustSpeech and Sig. Proc. 33(2) : 387-392, April 
1985. 



4. Estimer les retards 



= argtnin J£,^(t) 



XmU^L a(-)"Qi(R)<^i/")) ^^^^ 



Pour ce faire, appliquer le processus iteratif decrit ci-dessous : 

- Poser d = 1 . Obtenir une estimation f , de t^, en minimisant j^^(ti). 

t 

- Poser d = 2. Obtenir une estimation de en minimisant Ji/ln.ral) 

en fonction de et en consen^ant t, §gal a I'estim6e obtenue a la 
premiere etape. Puis, obtenir une nouvelle estimation de [to,.ToJ en 
mettant en oeuvre une minimisation multidimensionnelle de J£,<i(ln,r2l) 
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en utilisant les valeurs [fpfj] que Ton vient d'obtenir comme valeurs 
initiales. 

- Iterer ce processus jusqu'a ce que d = J, pour tous les chemins 
possibles de maniere a obtenir le vecteur des retards. 

5. Estimer la valeur du vecteur d'antenne par 

h = vp„,„(G.^.j(r)''Qi(R)Gi /x)). 

6. Former Testimee de la reponse Impulsionnelie : h = j(-r)a 

Precede sans a priori 

La deuxieme variante du prbcede correspond a un algorithme 
d'estimation parametrique totalement ou en majority autodidacte. 
L'algorithme associe tient compte, comme dans la premiere variante 
de realisation du caractere spdculaire du canal de propagation afin de . 
determiner ses parametres, notamment le vecteur d'antenne et le 
vecteur temporel caracterlstiques. Contrairement a la premiere 
variante, il he necessite pas la connaissance a priori du filtre 
d'emisslon. II suppose simplement que le filtre d'emission est a bande 
limitee et que sa bande est approximativement connue. 
Les differentes etapes du precede ou de Talgorithme sent par exemple 
les suivantes : 

1 . Estimer Lg la longueur du filtre d'emission.et clioisir L > -f Ar^ax et/c 

2. Appliquer la mdthode sous-espace decrite au point 2 du precede 
avec a priori : 

- Estimer ia matrice de covarlance R par exemple en mettant en 
25 geuvre les etapes decrites precedemment correspondant a 

l'equation(Eq. 28). 

- Galculer la matrice de projection sur I'espace bruit a partir des 
vecteurs propres associes aux pq(K + 1)-(K+L + 1) valeurs propres 
nulles de la niatrice, R . 

30 - Former la matrice Q£^(R) = Vi{nf^Vi^{fii)". 
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- Obtenir £[ le vecteur propr© assocle a la plus petite valeur propre 
d© la matrice Qx,(R)> 

et former %(^)- 

3. Former le nouveau critere param6trique : 

- Choisir la valeur de la longueur du canal k > i-x = pCL + 1) - i 

- et former la matrice 7>(hx). 

' Calculer la matrice contenant les vecteurs propres associes 

auxq(R+1)- " («+ ^i + i) 

plus petites valeurs propre de " Th(&i)"- 

- Former la matrice XJi COij.)- 

4. Estimer le nombre de trajets d , par exemple a I'aide de ia methode 
pr6c6demment referencee pour I'etape 3 de la premiere variante de 
realisation. 

5. Determiner une valeur minimal© de la bande du filtre d'emisslon : p, 
generalement p = 1/T. 

6. Choisir v (t), un filtre continu a bande limitee de bande p. Former le 
filtre v(t), 

v(t) = vit) 0<t<Lv<Lg ^^2) 
v{t) ^ 0 ailleuis 



7. Estimer les retards 

•f- = argi^n f^^(T) 

avec 
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ert appliquant le processus iteratif propose pour la methode sous- 
espace parametrique. (La matrice V^/r) est definle de la m§me 
fa9on que la matrice Gi/r) mais k partir du filtre v(t).) 
8. Estimer les vecteurs d'antennes : 

Applications possibles 

Le precede selon linvention trouve son application, par 
exemple dans les domaines mentionnes ci-apres, donnes a litre 
illustratif et nullement limitatifs. 
Rdseaux cellullaires 

Dans le contexle des transmissions cellulalres, la norme de 
transmission est connue. Le filtre d'emission normalise est done connu 
et il est pref6rable d'appliquer le proc6d6 avec a priori, utilisant la 
connaissance a priori du filtre d'6mission/r6ceptlon. Les applications 
potentielles sont les suivantes : 
Egalisation 

La propagation joue le r6le d'une convolution temporelle. A 
I'exception des communications a tres bas-d6bit (quelques centaines 
de bits/seconde), cette dispersion temporelle induite par la 
propagation cree de I'interference entre symboles (lES). Afin de 
pouvoir recuperer au mieux les symboles transmis, I'effet dispersif du 
canal doit etre compens§ par une etape d'egalisation telle que decrile 
dans I'ouvrage ayant pour titre "Digital Communications", pour auteurs 
J.G.ProakIs, publie aux Editions McGraw Hill, 1989. Si le filtre 
d'emission/reception est connu. k partir de la connaissance des 
parametres du canal de propagation, il est possible de reconstruire la 
r^ponse globale puis de mettre en oeuvre une §tape d'egalisation. 
L'int§ret d'une telle approche est qu'en general en estimant les 
parametres de la propagation plutot que la reponse globale. le nombre 
de parametres a estimer est r^duit, ce qui permet en general de 
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reduire la variance des estimateur. On peut ainsi amehorer 
significativement la qualite de lldentiflcation et les performances de 
I'egalisation. De plus, certains parametres caract§rlstlques du canal de 
propagation (les retards et les angles d'arrivee) 6voluent peu d'une 
trame a I'autre. La modelisation parametrique peut permettre de suivre 
de fagon plus efficace les variations des parametres du canal et 
d'am6llorer encore la qualite de fidentification. par exemple tel que le 
decrit la reference "Multi-channel MLSE equalizer with parametric FIR 
channel identification" de J.T. Chen et A. Paulraj. In Proc. VTC. pages 

710-713, 1997. 
Localisation 

La localisation de I'imetteur, problfeme classique des 
applications militalres. est devenue r§cemment d'un grand interet pour 
des applications civlles en communications mobiles. En effet, la 
localisation des mobiles dans un reseau cellulaire pemiet de rendre 
de nombreux sen/ices tels que la localisation des appels d'urgence. la 
gestion de la taxation, la detection des fraudes. I'aide k la planification 
cellulaire. I'amelioration de la mise en oeuvre du « handover », telle 
que decrite dans la reference » Position location using wireless 
communications on highways of the future de T.S.Rappaport. 
J H.Reed et B.D.Woemer, IEEE Commun. Mag., October 1996. 

II existe plusieurs techniques de localisation. La plus classique 
est la triangularisation. Dans ce cas. la localisation est e.ectuee k 
I'aide de deux ou trois recepteurs (stations de base) distincts estimant 
chacun les angles ou/et les temps d'amvee du signal emis telle que 
celle decrite dans la r§f6rence "Mobile Localization in a GSM Network 
"de H.El Nahas. PhD thesis ENST. june 1999. Ce systeme presente 
certains inconvenients qui ont pousse au developpement 
d'algorithmes de localisation a partir d'une station unique. Ces 
techniques de localisation n6cessitent toujours la connaissance a 
priori des caracteristiques de la propagation. Un precede de 



1 er depot 



19 

localisation est propose par exemple dans la reference " Precedes de 
localisation spatio-temporel de radlomobiles en rtillieu urbain " de M. 
Chenu-Toumier, A. Fen-^ol et J-J. Monot, Technical report, TCC, 
1997. 

Rltrage d'antenne - SOMA 

La connalssance des caract6ristiques spatlaies des differents 
trajets peut etre utilisee afin de focaliser remission dans la direction 
principale d'arrivee et ainsi minimiser I'effet des trajets multiples au 
niveau des recepteurs (uplink beamfomner design, downlink 
transmission design). 
Surveillance du spectre 

Dans le contexte de I'analyse radio-electrique du spectre, le 
rdcepteur n'a aucune information a priori sur I'emetteur. Le systeme de ^ 
reception scrute le spectre, Intercepte un signal transmis par un 
emetteur inconnu, et I'analyse afin d'en extraire certaines infomiations 
qui pemiettront d'ecouter remission ou de localiser I'emetteur. Dans 
ce contexte, il eist peu probable que I'on connaisse le filtre d'emission, 
et le precede sans a priori pour estimer les parametres du canal de 
propagation (. ; a) est plus adapte. Dans ce cas, le filtre d'emission 
n'etant pas connu, il n'est pas possible d'egaliser le signal regu, mais 
la connalssance des parametres (. ; a) permet par exemple la mise en 
ouvre des applications donnees ci-dessous. 
Localisation 

La principale application est la localisation. En HF, la 
localisation h station unique utilise classiquement une etape initiate de 
• provisions ionospheriques. Connaissant Tangle d'incidence et le profil 
de la densite 6lectronique (par provision ionospherique), une methode 
de trace de rayons utillsant la loi de Descartes pemiet de reconstituer 
le trajet de I'onde regue sur le rOseau de capteurs selon I'une des 
references suivantes : 
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. -Comparison of the .xing accuracy of single-station location of long 
range transmitters, de H.C.H'ring , IEEE Proceedings. 37(3) : 173- 
176, June 1990. 

. -loniospheric modelling in support of single station location of long 
range transmitters de L.F.MoNamara, Journal of Atmospheric and 
terrestrial Physics. 50(9) : 781-795. 1988. 

En outre, la connaissance des temps d'amvee pem^et de lever 
certalnes ambiguites dues k des multibonds. 

Autres . 
U connaissance des veeteurs d-antennes perniet de faire du 

flRrage spatial. U connaissance des retards pem,et de recombiner les 
traiets de propagation en phase. L'ullllsation dune des deux 
connaissances ou leur uMlisalion coniclnte peut constltuer une 
premiere «ape 4 la mlse en ceuvre d'un systSme d'enalyse du signal. 
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REVENDICATIONS 

1 '. Proced§ autodidacte ou partiellement autodldacte pour determiner 
les parametres spatio-temporels caracteristiques d-un canal de 
propagation dans un systeme comportant au moins un capteur de 
reception recevant un signal y(t) caract^rise en ce qu'il comporte au moins 
une etape oii la stmcture de type sp§culaira du canal est prise en compte 
et une etape de d§temnination conjointe des paranrietres, tels que les 
vecteurs d'antennes (a) eVou les vecteurs temporels (x) a partir des 
statlstiques d'ordre 2 des signaux regus. 

2. Proced6 selon la revendlcation 1 caracterise en ce qu'il comporte 
une 6tape de surechantillonnage du signal regu. 

3. Proc§d6 selon la revendlcation 1 caract6ris6 en ce qu'il comporte 
une 6tape ou les signaux sont regus sur au moins deux capteurs et une 
§tape de surechantillonnage du signal regu. 

4. Procede selon I'une des revendications 2 et 3 caracterisd en ce que 
la periode d'echantillonnage correspond a T/p avec T la p6riode symbole. 

5. Procede selon I'une des revendications 1 a 4 caracterise en ce qu'il 
comporte au moins les etapes suivantes : 

• choisir la longueur du canal global de propagation L>Lg + At^^ 

et la valeur du nombre d'obsen/ation K^L Avec Aw la plus grande 
valeur possible du retard relatif entre deux trajets, 

• determiner la matrice Q£(R) du critere sous-espace^ 

• estimer le nombre de trajets 4 
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estimer les retar(j$ . 

estimer la valeur du vecteur d'antenne iaar . ' ' ^» = vp„.i„(G£.j(r)^Qi(R)Gi 
former Testimee de la r^ponse jmpulslonne|l# : *^ = 

6. Proc§d§ selon I'une des rovendications 1 a 4 caracterise en ce 
u'il comporte au moins les ^tap0s suivantes : 
Estimer Lg la longueur du filtr© d'^mlssion, et cholsir 

L^Lg + AT^ et K^t, 

Appliquer la methods sous-espace d^crite au point 2 du precede avec 
apriorl: 

- Estimer la matrice de Qoyariance ^ . 

- Calculer la matrice de projection sur I'espace bruit n^. 

a partir des vecteurs propres assooies aux pq(K+1) - (K + 
vateurs propres nulles de la matrice R . 

- Former la matrice - Q^(R) = i?i(n£)Pt(«t)"- 

- Obtenir le vecteur propre associ§ h la plus petite valeur 
propre de la matrioe Qi(»-) , etfomrief hxC^) 

• Fonfner le nouveau crit^re param6trique : 

- Choisir la valeur de la longueur du canal R>Lt=p{.l + i)-i 
et fomner la matrice TM^^e) 

- Calculer la matrice tt^^ cbntensint les vecteurs propres 

p ■ ■ 



associes aux g(R + !)-(«+ + 1) 

P 
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plus petltes valeurs propres de 1> 
Former la matrice ViA^'ir) 

Estimer le nombre de trajets d , 

Determiner une valeur minimale de la bande du filtre d'6mlsslon : a 
Choisjr v(0 , un filtre continu a bande limitee sie bande b . Former le 

filtrev(t), iv{t)=Ht) 0<t<Lv<Lg 

\ v{t) = 0 ailleurs 

• Estimer les retards 

• Estimer les vecteurs d'antennes :. 

■fr p "» " 

7. Application du procede selon I'une des revendicatlons 1 a 6 pour 
la surveillance du spectre d'un canal de propagation k des fins de 
localisation a partir d'une ou de plusieurs stations en HF. 

8. Application du procede selon Vune des revendicatlons 1 h 6 dans 
le cadre de liaisons de communication normalisees k r§galisation ou la 
localisation ou encore le filtrage spatial. 
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